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I. INTRODUCTION 

Light-scattering has proven to be an important tool 
for investigating condensed matter physics. In the field 
of supercooled liquids, the structural relaxation covers 
many decades either in the time, or in the frequency 
domain, the latter being accessible by e.g., Fabry-Perot 
techniques. The measured spectra, see e.g. reflect the 
slowing down of the structural relaxation upon lowering 
the temperature and exhibit the nontrivial power-laws 
and stretching effects found by other techniques, such 
as dielectric spectroscopy for instance. The most direct 
measure of the hindered motion due to the cage effect can 
be observed by depolarised light-scattering in the back- 
scattering geometry. For depolarised light-scattering per- 
formed at other scattering angles, one observes an admix- 
ture of the transverse current motion to the pure back- 
scattering signal. Similarly, for polarised scattering, one 
obtains a contribution from density fluctuations which 
result in the Brillouin resonance ||^. 

The subtlety of light-scattering lies in disentangling 
the dependence on frequency shift, oj, and wave- vector 
transfer, q, as well as on incident and outgoing polari- 
sations. Since the wave- vector transfer, q, is small for 
light-scattering, a generalised hydrodynamics approach 
is suitable. There, the spectra are described in terms of 
a number of frequency-dependent memory kernels, e.g. 
viscosities. These kernels have sometimes been written 
on the basis of heuristic arguments. This is the case, 
for instance of in which different previous attempts 
are also described and discussed. A more fundamental 
approach consists in deriving them from a microscopic 
theory through, say, a Zwanzig-Mori technique. The 
first such attempt was made by Andersen and Pecora , 
who, in fact, made a purely formal use of the technique, 
the memory kernels being eventually approximated by in- 
stantaneous interactions (Markov approximation). Much 
more recently, the technique was used in its full generality 
in 1^. Using only general symmetry considerations, |^ 
showed that the description of the light scattering spectra 
involved ten frequency-dependent functions. This large 
number was the price to be paid in order not to miss any 



effect that leads to a fluctuation of the dielectric tensor 
(5e.y(q, t). In the present paper, we shall follow an in- 
termediate route; we shall derive, for a selected set of 
dynamical variables, the precise form of their equations 
of motion, and of the corresponding relaxation kernels. 
Our results are valid whatever the temperature but are 
restricted to the case of molecular supercooled liquids 
formed of symmetric top molecules. Their application 
to the light scattering problem requires a precise form of 
Seij{r,t); following Eq. (12) of Part I of this series of 
papers, we shall assume 5eij{q,t) to depend only on two 
variables of the problem, namely the density and the ori- 
entation fluctuations. Then the light-scattering problem 
is reduced to calculating the density-density, orientation- 
orientation as well as mixed correlation functions which 
are expressed with the help of appropriate memory ker- 
nels. 

As a further simpliflcation, we shall ignore temperature 
fluctuations, i.e. the hydrodynamic poles associated with 
energy conservation. This restriction is probably justified 
for Brillouin scattering experiments, since, for the scat- 
tering vectors involved, the Rayleigh line lies at so low 
a frequency that it is inaccessible to the usual frequency 
domain methods. Furthermore, for liquids, the ratio of 
the isobaric heat capacity to the isochoric one is close to 
unity and, correspondingly, the total weight of the Bril- 
louin lines is much larger than the weight contained in 
the Rayleigh line; similarly, the isothermal sound veloc- 
ity is close to the adiabatic one. The situation is different 
for time-based methods like impulsive thermal stimulated 
Brillouin scattering 0, |^ where the heat diffusion con- 
tribution can be observed as a late stage of the relaxation 
signal, but this aspect of the problem will not be dealt 
with here. 

The goal of this paper is twofold. The first is to give 
a microscopic derivation of the constitutive equations for 
the density and orientation fluctuations used in |^, 

, and to derive some new results from this microscopic 
approach The second is to compare the results one 
can obtain from the three microscopic approaches |^ 
and the present paper, which differ in the variables taken 
into account and/or in the scattering model. 

Consequently, this paper is organized as follows. The 
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phenomenological equations of ^ are microscopically 
derived in Section II. In particular, we show that the 
four memory functions which enter into those equations, 
namely the bulk viscosity, rjt (t) , the center-of-mass shear 
viscosity, r]s{t), the rotational friction T'{t), and the 
rotation-translation coupling, can be expressed in 

terms of the dynamical variables of the problem, and 
of a reduced time evolution R'{t) which does not con- 
tain the hydrodynamics poles of the problem. Similarly, 
the rotation-translation coupling constant. A', and the 
molecular libration frequency, Uq, which are the other in- 
gredients of these equations of motion, will be expressed 
in terms of equal time thermal averages of some vari- 
ables of the problem. We make use of these microscopic 
expressions of the memory functions in Section III to de- 
rive necessary conditions on the imaginary part of their 
Laplace transform, and on some contributions of them. 
These conditions will be such that the light scattering 
spectra will be always positive whatever the values of 
the coefficients linearly coupling the density and orien- 
tation fluctuations to Je^ (q, t) . Section IV makes use 
of the same expressions of the memory functions to re- 
late, through a Green-Kubo formalism, the correlation 
functions of some variables to specific combinations of 
the Laplace transforms of these memory functions. In 
particular, we shall show that rjtit) can be expressed as 
such a correlation function. The same will be true for 
r]T{t), the memory function which takes into account all 
the retardation effects related to the propagation of the 
transverse phonons; this is not a priori obvious because 
r]T{w) will turn out to depend in a complex way on the 
Laplace transforms of several memory functions defined 
above, as well as on A' and luq. Section V will compare 
the expressions for the light scattering intensities that can 
be obtained using the three sets of variables and of di- 
electric fluctuation models already mentioned. We shall 
show that the set proposed in Q leads to awkward forms 
of the relaxation kernels when they are not restricted to 
a Markov approximation, but used for a molecular su- 
percooled liquid. Conversely, as expected, the results 
obtained in Part I are a restriction of those of j|] cor- 
responding to definite simplifying assumptions. A brief 
summary and some comments conclude the paper. 



II. A ZWANZIG-MORI DERIVATION OF THE 
DYNAMICAL EQUATIONS 

We consider a dense liquid of N linear molecules of 
mass m at temperature T enclosed in a volume V. Sta- 
tistical correlations of phase space variables in terms of 
the Kubo scalar product futl, {A{t)\B) ^ {SA{t)*SB), 
SA = A — (A) , provide the simplest information on the 
system's dynamics with (.) denoting canonical averaging. 
The thermodynamic limit, N —^ oo, with fixed particle 
density, n = N /V , is implied throughout. The time evo- 
lution of the observables is driven by the Liouvillian C: 
dtA — iCA = {H^A], where H denotes the Hamilton 



function and { , } the Poisson bracket. We consider the 
dynamics of the fluctuating molecular orientation tensor, 
written directly in the reciprocal space: 

g,, (q) = [uo^.u^, - j e^i-^° , (1) 

where the degrees of freedom of the a-th molecule are 
specified by a unit vector, u^, for the orientation and by 
the position of its center-of-mass, R^. The spatial modu- 
lation of a fluctuation is characterized by its wave vector, 
q, and latin indices denote cartesian components. The 
9 components of Qy (q) are not independent, since the 
orientation tensor is symmetric and traceless, reducing 
the number of independent components to 5. The nor- 
malisation is chosen such that the correlation functions 
are intensive. 

Furthermore, we consider the fluctuations in the mass 
density: 

N 

p(q) =miV-i/2^exp(iq-R„) , (2) 
and the cartesian components of the mass current: 

N 

J,(q) = TV- ^ exp (iq • R„) , (3) 

a=l 

where denotes the momentum of the a-th molecule. 
(Equivalently, one could use the particle density n(q) = 
p(q)/m and the velocity Wi(q) = Jiiq)/pm, where pm = 
mn is the mean mass density). 

A. Static averages 

The static correlation functions need to be evaluated 
to lowest order in q only. Since the Hamilton function 
respects rotational invariance, all static averages in the 
liquid phase have to remain unchanged under any rota- 
tion of the system: this implies that, e.g. correlators 
between, say, any second rank traceless tensor and any 
scalar will vanish in the long-wavelength limit (see, e.g. 
Eq. (0)). 

The static average of the density can be expressed as: 

(p(q)|p(q))=mV/c2 + 0(g2), (4) 

where c is, here, the isothermal sound velocity |]l8|| de- 
fined in terms of the long-wavelength limit of the static 
structure factor via = v'^/S{q — > 0), while v = 
^Jkjfrjra denotes the thermal velocity. As usual, the 
current correlations read: 

(J,(q)|Jfe(q)) = (5,fcmV. (5) 

To lowest order in g, the equal-time correlators of the 
tensor variables read: 

(Qz,(q)|Qfci(q)) = ^'A„-,M + 0{q^) , (6a) 
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where: 



^ij,ki = ( SikSji + SilSjk — -SijSki 



(6b) 



is a fourth-rank tensor, the structure of which is governed 
by rotational symmetry. The long-wavelength limit of 
the 9x9 correlators in Eq. ( |6a]) is thus determined by a 
single number, S'^, denoting the long- wavelength limit of 
the corresponding generalised structure factor, a quantity 
which is, as in Eq. (^), proportional in leading order to 
keT. 

Due to rotational symmetry, the overlap of the ten- 
sor variables with the density vanishes in the long- 
wavelength limit according to: 



(g.,(q)|p(q))-0(g2 



(7) 



We shall also need to consider the tensor currents, 
Qij (q) = iCQij (q) , which are normalised by: 



(Q^j(q)IQfci(q)) -^^'A 



ij^kl 



(8) 



with the characteri stic frequency scale f2; we show, in 
Appendix A, Eq. (^), that fl = y^2k^Tj5I, where / 



is the moment of inertia of the molecule for a rotation 
around an axis perpendicular to the molecule symme- 
try axis and passing through its center of mass. The 
ratio of the static averages of the orientation and the 
orientational current, loq = Q/S, will determine the ax- 
ial libration frequency, a frequency characteristic of the 
short-time expansion for the orientation correlation func- 
tion (see Eq. (27c)). Hence, there is a close analogy be- 
tween the set of Eqs. ( |6a[ ) and (^) and the set of density 
plus momentum current correlators whose ratio deter- 
mines the isothermal sound velocity c characteristic of 
the initial decay of the density correlators; c and ujq are, 
to leading order, independent of temperature. 

The correlation function between the mass current and 
the tensor current components has now to be considered; 
it is strictly equal to zero, whatever is q: 



(O.,(q)|Jfe(q)) = 0. 



(9) 



This is due to the fact that we put the point of refer- 
ence of each molecule, Rq,, at its center-of-mass, (see 
Appendix A for a thorough discussion of this property). 
The remaining static averages between the four distin- 
guished variables p{q) , Ji (q) , Qij (q) and Qij (q) vanish 
due to time reversal symmetry. 



B. Constitutive Equations 

The mass conservation law relates the density to the 
momentum current: 

dtp{ci,t) = iqkJk{<\,t) . (10) 

Similarly, the conservation of momentum yields: 

dtJk{ci.t) ^ iqiYVuici.t) , (11) 



where 11^; (q, t) denotes the fluctuating momentum cur- 
rent tensor. At last, we can write the trivial identity: 



afgjj(q,i) = (3,j(q,t), 



(12) 



which defines Qij{(\,t) as an orientational tensor force. 
In order to close the system, we need constitutive equa- 
tions for the momentum current tensor, nfc;(q, t) and the 
orientational tensor force. This will be achieved here 
through generalised Langevin equations which will intro- 
duce appropriate memory kernels. Let us first introduce 
the projection operator, P: 



P = 



l<3fci(q))^(<3fei(q) 



-\Qkim^{Qki{ci)\ + \p{ci)) 



-kfc(q)) 



1 
1 



(p(q)l 



(Jfc(q)|+0(g2), 



(13) 



where the sum over repeated indices is implied. P is a 
projection operator because, once the symmetric charac- 
ter of Qy (q) and (3ij(q) has been taken into account, 
one can check that indeed P^ = P. P projects onto 
the subspace spanned by density, mass current and the 
symmetric traceless parts of the orientation and the cor- 
responding current. 

The time evolution operator, R{t) = cxp{iCt), can be 
exactly reformulated as 



R{t) = R{t)P + [ R{s)PiCR'{t - s)ds 
Jo 



+ R'{t), (14) 



with the reduced operator R'{t) = Q exp{iQCQt)Q, 
where Q = 1 — P, and a short proof of Eq. ( p^ ) is given 
in Appendix B. The benefit of this procedure lies in the 
following. The time evolution operator, R{t), possesses, 
in addition to a non-hydrodynamic part, long-lived hy- 
drodynamic modes that are due to conservation laws. 
This leads to resonances in the spectra, viz. the Fourier 
transforms of the time correlation functions of all the dis- 
tinguished variables for small but nonzero wave vectors. 
Conversely, the reduced time evolution operator, R'{t), 
devoids the hydrodynamic singularities and correlation 
functions with R'{t) are regular in the long- wavelength 
limit. The problem of handling the slow relaxation due to 
hydrodynamic conservation laws is treated explicitly in 
the low-dimensional subspace of the distinguished vari- 
ables. On the contrary, the slow structural relaxation will 
be dealt with the help of correlation fimctions of R'{t), 
the second term of the r.h.s. of Eq. (|T^), which will 
appear in the form of memory kernels. In the spirit of 
generalised hydrodynamics, the long-wavelength proper- 
ties are described properly by keeping the wave- vector de- 
pendences introduced explicitly by the conservation laws, 
while the memory kernels can be evaluated in their long- 
wavelength limit. 

Before deriving the constitutive equations from Eq. 
(O) for the missing quantities, 11^ (q, i) and Qij{ci,t), 
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some comments on the structure of the resulting equa- 
tions in the long- wavelength limit are in order. First, 
from time reversal symmetry, the instantaneous coupling 
(first term of the r.h.s. of Eq. (|l^) will be non-zero 
only for variables of identical time-parity. Since both the 
momentum current and the oricntational force have even 
time parity, this instantaneous part will consist of den- 
sity and orientation only. Second, rotational symmetry 
implies that the coupling of irreducible tensors of differ- 
ent ranks is suppressed in the long-wavelength limit by 
appropriate powers of the wave number. The dynamical 
correlators enjoy the same property, since the time evo- 
lution does not change the rank of a tensor. We shall 
keep only the lowest nontrivial terms in this small-wave- 
number expansion, as was already hinted at by keeping 
only the lowest order of the static averages in the pre- 
ceeding section. 

Finally, in the second term of the r.h.s. of Eq. (14), 



one can let iC operate on the 'bra' part of the projector, 
P; for instance: 

|p(q))^T(/5(q)l*^ = -|p(q))^rT(p(q)l 



(15a) 



because R'{t — s) contains, on its l.h.s, a Q = 1 — P factor, 
the contribution of Jfe(q), and similarly of Qki{(l), are 
eliminated from this second term and one obtains: 



1 



PiCR'it - s) = -|Q^i(q))_(Qfc;(q)|P'(t - s) 



mKB-L 



(15b) 



Let's first handle the momentum current tensor which 
we decompose into: 



% (q, t) = 6^jp{q, t) + Tiij (q, t) 



(15c) 



Here: 



p(q, t) - [n,,(q, t) + T\yy{ci, t) + n,,(q, t)]/3 (15d) 

denotes the fiuctuating pressure so that 7ry(q, i) is a 
traceless symmetric second rank tensor. Multiplying Eq. 
( p^ ) from the right by p(q) yields the desired Langevin 
equation for the pressure fluctuation: 

p(q,i) = R{t)Pp{q)+ [ R{s)PiCR' {t ~ s)p{q)ds 

+R'it)piq) . (16) 

The first term represents an instantaneous coupling to 
the distinguished variables of the projector. The second 
yields a retarded coupling and the last term is a rapidly 
fiuctuating term that we shall call 'noise', i.e. which is 
uncorrclated for all times to the distinguished variables. 
Hence, this term can be ignored for the evaluation of 
the correlation functions of the distinguished variables. 



Nevertheless, the same term will be useful in establishing 
the Kubo formulae of Section IV, which will relate the 
time-dependent correlation functions of some variables to 
the memory kernels of the dynamical equations. 

In order to evaluate the first term of Eq. ([l6|), we 
need static correlations of the pressure with the distin- 
guished variables. Time- inversion symmetry allows non- 
vanishing correlations only with the density and the ori- 
cntational tensor. Since rotational invariance implies: 



(7ri,(q)|p(q))=0(<z2), 



(17a) 



one can evaluate (p(a),p(c[)) by using the conservation 
of momentum, Eq. (ll|), and Eq. (|5|), up to terms of 
order 0{q'^): 

(p(q)b(q))«% = (/o(q)b(q))i9/4; 

= (p(q)N-Z/Hfe,(q)) = (p(q)|jfe(q)) 

= -(p(q)|Jfe(q))=z<Zi(Ji(q)|Jfe(q)) 

(17b) 



This yields: 



(p(q)|p(q)) = "^V-fO(q2). 



(17c) 



Conversely, rotational invariance implies, similarly to 
Eq. d), that (p(q)|Q,j(q)) = 0{q^). Collecting the 
terms appearing in P, the long-wavelength instantaneous 
coupling is then simply given by: 

P(t)Pp(q) =cV(q,t). (18) 

Let's turn now to the retarded couplings. From Eq. 
( |l5b| ), one has to consider only the couplings of p(q) 
with QfcKq) and with JA;(q), and these couplings involve 
{QuMWit - s)b(q)) and (Hfe,(q)|i?'(t - s)|p(q)). As 
Qfc;(q) is traceless while Iikl{^) is not, the first term is 
0{q^) while the second is of order unity. Consequently, 
to lowest order in g, we need to consider only the re- 
tarded coupling to Jfc(q). It is convenient to introduce 
the long- wavelength pressure correlator in the form: 



{p\R'{t)\p) 



= Vbit) . 



(19) 



(Here, and in the rest of the paper, omitting the wave- 
number dependence indicates that the quantity is to be 
evaluated for q ^ 0). 
One thus obtains: 



iqiiUki{<D\R'it^s)\piq)) 
= iqkVb{t- s)kBT/n + 0{q^) 



(20a) 



once we make use of the fact that (TTfc; (q)|P'(t)|p(q)) is 
of order 0{q^). Introducing Eq. (20a) into the expres- 
sion of the retarded coupling and replacing, in Eq. (p^), 
Jki<i) by mnvk{(l), one obtains the standard, retarded, 
constitutive equation for the fiuctuating pressure: 

p{q_,t) = c^p{c[,t)+i / rib{t ~ s)qkVk{q, s)ds 
Jo 

+noise , (20b) 
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which is the Fourier transform of the usual equation: 



Similarly: 



p(r,t) 



Jo 

-\-noise . 



s) div v(r, s)ds 

(20c) 

Let us now apply the same technique to 7rij(q), the 
traceless part of Iltj (q) . Multiplying Eq. (|lj) iDy tt^ (q) 
from the right yields the Langevin equation for the trace- 
less part of the momentum current tensor. In order to 
evaluate the instantaneous couplings, one needs to com- 
pute the static correlations of this tens or w ith the dis- 
tinguished variables. Because of Eq. (17a) and of the 



time reversal symmetry, one is left with the sole evalu- 
ation of (Qfc;(q)|7rij(q)) for which neither tensorial nor 
time reversal symmetry constraints apply, in the long- 
wavelength limit. However from momentum conserva- 
tion: 



«9i((3fci(q)kjj(q)) = (Qfei(q)|i/:j.,(q)) 
= -(Qfci(q)l^z(q)), 



(21a) 



where the first equality is correct up to order 0{q^) due 
t o th e traceless character of Qfc;(q). The r.h.s. of Eq. 
(21a) is equal to zero at every order in q, Eq. (^, so 
that there is no instantaneous coupling of 7rij(q) with 
the distinguished variables, in the leading order in q con- 
sidered in the present paper. Hence, one is left with the 
evaluation of the memory kernel, which splits into two 
parts: 

-the traceless momentum current tensor auto- 
correlator which, in line with Eq. (|l9|), we write in the 
form: 

71 

i7rki\R'{t)\7T,j) ^ = Vsit)\j,ki , (22) 

and we call rjsit) the time-dependent shear viscosity, as 
it couples the momentum current to the strain rate, see 
Eq. (||). 

-the coupling of the traceless momentum current with 
the corresponding orientational force, that we write as: 

'Qki\R'{t)\TT,j) = -fi{t)A,jM ■ 



(23) 



For reasons similar to those used in Eq. (20a) 

iqiiUui(i)\R'{t~s)\n,,{q)) 
= iqi{nkiiq)\R'{t - s)|7ry(q)) 



iqiVsit - s 



^kl,ij 



(24a) 



(24b) 



and: 



iqi'^^^^l^ki,ij = im[qiVj{(\} + qjVi{q) - ^(5ygfeVfe(q)] 



= -rriTij (q) , 
where Tij(q) is the strain rate tensor, so that: 
1 



(24c) 



i?(s)|Jfe(q)) 



mkBT 



^g^(^fc^(q)|i?'(^-s)|x,,(q)) 



'Ris)\Qkli^))^{Qkl{cl)\R'{t ~ s)K.,(q)) 
R{s)\Qki{(l))^Aki,;jfi{t - s) 



Qy(q, s)^i{t - s) . 



(24e) 



Then, the generalised constitutive equation for 7r,;j (q, t) 
reads: 

7rij(q, i) = - / ris{t- s)T^j{q,s)ds 
Jo 



+ / fj,{t — s)Qij{q, s)ds + noise . (25) 
Jo 



Combining Eq. ( pOb| ) and Eq. (pSj), one obtains: 

njj(q, i) = SijC^p{c[,t) +iSij / T]t{t - s)qkVk{q, s)ds 

Jo 



ris{t~ s)Tij{q,s)ds 



+ / fj,{t — s)Qij{q, s)ds + noise , (26) 
Jo 

which is exactly the Fourier transform of Eq. (3) of Part 
I, once one has noted that H^ (q, t) is the opposite of the 
stress tensor, aij{q,t). 

To derive an equation of motion for Qij{q), we again 
make use of Eq. (|l^). For the instantaneous contribu- 
tion, R{t)PQij{q), only the term involving the orienta- 
tion in the projector P needs to be considered, the other 
terms dropping out for tensorial or time-reversal symme- 
try considerations. Since: 

{Qkl{^)\Q^J{^)) = -(QkiM^.iq)) 

= -^^2Afe^,^J, (27a) 

i?(OPQy (q) = -UJ^„Q^J (q, t) , (27b) 

with the axial libration frequency: 

LUo = n/S. (27c) 

The evaluation of the retarded couplings proceeds along 
the same lines as for the momentum current tensor. 
Defining: 



(QkllR'itmj)-^ - T'{t)Akl,^,, 



and: 



A' 



n'^n _ 2n 
~ 51 



(28a) 



(28b) 



j{q,s)r],{t- s). 



(24d) 



one easily obtains, with the help of Eqs. ( |l5b[ ), and (24c): 

PiCR'it-s)Q,j{q) = -|Qy(q))r'(i-s) 

+A'|r„(q))/i(t-s), (28c) 
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once one has noted that, because of the traceless charac- 
ter of Qij{q): 



{nki\R'{t-s)\Q^,) = {Trki\R'{t-s)\%j). 
Collectmg the various terms, one thus obtains: 



(28d) 



(q,i) = -ujlQ^j{q,t) - / T'{t- s)Q,j{(i,s)ds 



A' / n{t - s)nj{q,s)ds 

JQ 



(29) 



The same memory kernel, fi(t), occurs in the constitu- 
tive equation for the orientational force, Eq. (|2g|), as 
response to a momentum gradient, and in the equation 
for the momentum current, Eq. (25), as a reaction to an 
orientational current. This can be considered as a general 
consequence of Onsager's principle, and it appears, here, 
naturally as the result of the use of the Zwanzig-Mori 
formalism. Equation ( p9| ) is, as expected, the Fourier 
transform of Eq. (4) of Part I, as briefly argued in [p^ . 
The Zwanzig-Mori formalism thus leads to the micro- 
scopic derivation of the equations proposed, on a phe- 
nomenological basis, in ^. There are, nevertheless, 
alrea dy tw o bonuses. One is the precise definitions of A', 
Eq. ( p8h| ), in terms of quantities a priori known, and 
of Wo, Eq. (27c), which can be obtained from thermal 



averages of {Qij\Qki) and {Qij\Qki)- The second bonus 
is the precise definitions, through R'{t) of the four mem- 
ory kernels, T]b{t),rj,{t), fi{t) and r'{t), Eqs. (|l|,||j23|) 
and (28a). We shall show, in the next two Sections, that 
these expressions allow: 

- on the one hand (Section III) to precisely define under 
which conditions, all the Brillouin intensities, derived or 
recalled in Part I, are positive whatever the frequency, 
within the scattering model used in Q. 

- on the other hand (Section IV) to show, through 
Kubo's formulae, that these kernels can, directly or indi- 
rectly depending on which one is considered, be measured 
as correlation functions of q ^ dynamical variables. 



III. THE ONSAGER RELATIONS AND THE 
POSITIVENESS OF THE SPECTRA 

A. Summary of the light scattering results of Part I 

In Part I, pl, making use of the equations of motions 
(Eqs. (p^, (pT|), ( p6| ) and (|29|)), we gave an expression 
for the intensity of the VV light-scattering spectrum un- 
der the assumption that the fluctuations of the dielectric 
tensor could be written as the linear combination: 



ey (q) = aSijp{q) + 6Qy(q) . 



(30) 



Taking the convention that the Laplace transform of f{t) 
would be /(w) = LT[f{t)]{Lj) = i dtf{t) exp{~iujt), 



this intensity was expressed (see Eq. 36, Part I) in 
terms of all the quantities defined in Section II, and of 
(|g"^,|2) = 5-2, see Eq. (Q. Using Eqs. (|27^ and 
( ^8b| ) which relate wo, A', il and S, the result obtained in 
Part I can be cast into the form [n9|: 



/yy(q, iv) — —Imi 

CO 



462 ^2 



3 D{uj) 



2hK' 
3mn 



1 2 



r{uj) 



(31a) 



where PL{q,Lu) is the longitudinal phonon propagator: 
PL{q,LLiy^ = uj^ - q^c^ - q^ujr]L{uj)/mn , (31b) 

with: 

4 A' 2 
7?l(w) = r]t{uj) + D{uj)r{uj) ] 



6 u) 
4 



D{lu) = ui'^ + ujT' {uj) - 



r(w) 



(31c) 



(32) 



(33) 



rixit), defined through the last line of Eq. ( pig ) is what 
we shall call the transverse viscosity. 

Similarly, the expression for the intensity of the VH 
light-scattering spectrum, already derived in jTl 0] within 
the same model, was recalled in Part I (Eq. (48)); with 
the present notations, it reads: 

Wq,c.) = -/m 



+q' 



A'r(cj) 



cos^ -m^v^PT{q,uj) \ , (34a) 



where 9 is the scattering angle and: 

Pt(9,w)^"'" — uP' — q'^ijjriT{tjj)/mn (34b) 

is the transverse phonon propagator. It is convenient to 
separate out the angular contribution in Eq. ( |34a| ) by 
rewriting it in the form: 



IvHi<l,uj) = 62 
with: 



sin2 -Ibd{^^) + cos^ -lTiq,uj) 



D{uj) 



(34c) 
(34d) 
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— Inil —— — r 
uj I D{uj) 



I 2 2 2 

+q m V 



'K'r{uj)' 


FT(g,a;)| 


mn 





(34e) 



B. Necessary conditions on the memory kernels 

In the present part of Section III, we show some general 
properties of the four memory kernels r]i,{t),r]s{t), ^{t) 
and r'{t) that can be derived from their microscopic ex- 
pressions. These properties are of interest for the light 
scattering spectra and, in particular, they are such that, 
when fulfilled, the spectra are positive whatever the fre- 
quency and the ratio b/a of Eq. ([3C|). 

Firstly, from their definitions, Eqs. (jl^, p3|, 28a), 



one checks that the memory kernels are real and have 
even time parity. Furthermore, by taking the special lin- 
ear combinations: 



Qo 



[20. 



Tryy]/Vl2, (35a) 
-Qyy]/y/l2, (35b) 



r]s(t) and T'(t) can be written as auto-correlation func- 
tions similar to r]{,(t): 



T\t) 



(^20|i?'(0k20), 



2(Qo|i?'(t)|Oo). 



1 



(36a) 
(36b) 



Thus the Laplace-transforms of the memory kernels have 
the usual properties described, e.g., in jll], ^ ^ (see 
also Appendix C). In particular, this Appendix shows 
that these Laplace transforms are analytic in the lower 
complex half plane and that the inequalities: 



Imrib{uj) > , 



Imr]s{uj) > , 



ImT'{uj) > 0, 

hold for all complex lu with ImuJ < 0. 
The mixed correlation function: 

~n^^l{t) = iQo\R'{t)\7T2o) 



(37a) 



(37b) 



(37c) 



(38a) 



can be read as an off-diagonal element of the matrix cor- 
relator built on Qo and 7^20- Since the imaginary part of 
the Laplace transform of thi s m atrix is positive semidef- 
inite (see Ap pend ix C, Eq. (|C6|)), one obtains, with the 
help of Eq. (|28b|), the inequality: 



[Im7js{uj)][Imr'{uj)] - A'[Imn{uj)]^ > 0, (38b) 



for all u! in the lower co mple x half- plane. The system 
of inequalities Eqs. (p7 



38b) is a generalisation of On- 
sager's relations to finite frequencies: from the micro- 
scopic approach, one obtains that the imaginary part of 
the matrix of kinetic coefhcients is positive definite for 
any frequency. 

Let us mention one useful consequence. First by 
Fourier back-transform: 



AV(t)2 = A' 



< A' 



< 



\Im^{uj)\ 



cLu 



ImT'icj) 



duj 



,(38c) 



where in the last line we used Eq. (38b). The last relation 
implies that the translation-rotation coupling is bounded 
by: 



A'fi{tf <rj,it = 0)T'it = 0). 



(38d) 



Positiveness of the light scattering spectra and 
further relations 



The four inequalities, Eqs. (37,38b), turn out to be 
suffic ient t o pr ove that the light-scattering spectra, Eqs. 
( |31a| ) and (34a), are positive for any frequency. The proof 
will be given for real frequencies only since the algebra 
greatly simplifies. By similar methods, one can extend 
the proof to hold for all frequencies in the lower com- 
plex half-plane. Let us first recall that if A is a sym- 
metric complex matrix, then one can write /to(A~^) — 
— [ImA+{ReA){ImA)~^{ReA)]^^. If ImA is a positive 
definite matrix, one proves that {ReA){ImA)^^ (ReA) 
has the same property, so that —Im{A~^) is also a pos- 
itive definite matrix; in particular, its diagonal elements 
are positive |2^. Let us now make use of this property 
to prove the positiveness of th e spe ctra, starting with the 
VH spectrum written as Eq. (34c). 

First: 



Ibd{<^) = 



ujllmT'{Lo) 



-ujReT'{uj)Y + [ujImT'{uj)]' 



(39a) 



so that Ibd{^) is always positive. Second, consider the 
matrix: 



LuD{uj) qD.ijj ijl(uj) / (mv) 

qVluj ii{uj) I Imv) —uj-\'q^rjs{uj)/mn 



(39b) 



(here the matrix elements we are not interested in are 
abbreviated by asterisks). One checks that, for real w, 
the imaginary part of the matrix on the r.h.s is positive 
definite: indeed, its diagonal elements are positive (Eqs. 
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37b, 37c) while the corresponding 2x2 determinant is 
proportional to the l.h.s. of Eq. ( |38b| ). Consequ ently , 
so docs the imaginary part of the l.h.s. of Eq. (39b). 
In particular, the diagonal elements on the left-hand side 
are positive. After some algebra one finds that: 



ImFT{q,uj) — Iriq,^) ' 



(39c) 



9 9 9 
TTl V 

ImCriq^uj) = ImPT{q,uj). (39d) 



As both Ibd{ijj) and lTiq,uj) are positive whate ver m 
real, the depolarised hght-scattering spectrum, Eq. ( |34c) ) 
i s alw ays positive. Also, from the form of PT{q,oj), Eq. 
(341:), the sign of ImPT{q,u!) is the same as that of 
77t(u;); this implies: 



Im rix{io) > . 



(39e) 



Thus, in spite of its intricate expression, Eq. ( 31c ), one 
can prove that ImrjTiuj) is always positive, a result which 
will be obtained again through the Green-Kubo tech- 
nique in Section IV. 

For the polarised spectrum, let us consider the similar 
matrix: 

FL{q,uj)/VL^ ujGLiq,Lo)/{nmvq) 
ujGL{q,i^)/{^mvq) uj'^CL{q,i^)/{'m^v'^q^) 



LoD{uj) 



(40a) 



Again, the imaginary part of the matrix on the r.h.s of 
Eq. (40a) is positive definite; its diagonal elements are 
positive, Eqs. (^7|), while the corresponding 2x2 deter- 
minant is proportional to: 

4 4 
[ImT'{uj)][Imr]b{uj) + -Imr^,{uj)] - -K'[Im . 

(40b) 

This ter m is also positive, because of Eqs. (38b) and 
( |37a , 37c). Explicit evaluation of the inverse of the matrix 
on the r.h.s. of Eq. (40a) yields: 



ImFLiq.uj) 



—Imi -, — - 

uj l3D(w) 



2 K'r{uj) 



3 mn 



m''v^PL{q,uj)Ym 



J r< I ^ 2 2 2 r f 2 AV(tj) ^ , , 

ImGL{q,uj) = q-'m^v^Imi —PL{q,Lo] 

u) 6 mn 



(40d) 



2 2 2 
Q TTl V 

ImCL{q,uj) = ^ ImPL{q,co) . (40e) 

to 



Si nce t he imaginary part of the matrix on the l.h.s. of 
Eq. (40a) is positive definite, this is also true for the ma- 
trix whose elements are ^^(q, w), Gl((J, w) and CL(g,u;). 
In consequence: 

2 

/l((j,cj) = a'^ImCL{q,oj) — -j=ab ImGLiqjio) 



Since 



V3 

1 

-ImFL{q,uo) > 0. 



ivviq,^) = b^iBoi'^) + h{q,^) , 



(40f) 



(40g 



the VV spectrum, Eq. ( 31a ), is positive, whatever a and 
b. [Note that the same technique could be applied to 
the HH spectrum, Eq. (43) of Part I, to prove that it is 
also positive, whatever the scattering angle.] Equation 
(40c) can also be used to prove that ImriLjoj ) > 0, a 
conc lusion which already resulted from Eqs. (B9q) and 
( |37a| ) . Let us stress that the positiveness of /yy (g, lu) for 
any uj, whatever q, is not a trivial result as Ihiq, w) is the 
sum of a g-independent term, proportional to Ibd{'^), 
and of a g-dependent term. One could naively think that 
Ivviq, '-^) could be always positive only if the same would 
be true for this g-dependent term. Figure 3 of Part I, |3[| 
shows that this is not the case. In fact, though Eq. (p8q) 
does not invoke g, it insures that ImT'{uj) is large enough, 
whatever to, for the sum of the two terms of lL{q,Lo) to 
be positive, independently of the value of g. A similar 
argument holds for the Ivniq,^) spectrum. 



IV. THE GREEN-KUBO APPROACH TO THE 
MEMORY KERNELS 

A. Preliminary Remarks 

The comparison between measured B rillou in sp ectra 
and their predicted intensities (Eqs. c 
requires the knowledge of the four memory 
riii{t),r]s{t), ii{t) andr'(t). Although microscopically well 
defined, those cannot be evaluated exactly, so that in 



and 34a) 
kernels 



practice, they are frequently taken as empirical fit func- 
tions. Whereas the correct memory kernels are guaran- 
teed to reflect all the restrictions of the correlated motion 
of the translational and orientational degrees of freedom, 
i.e. autom atically fulfill the relations, Eqs. ( |37a , 371: ,37c) 
and (38b), this needs not be true for these empirical func- 
tions. Hence, one has to carefully choose their parame- 
ters so that these relations are fulfilled. 

A possible intermediate approach consists in obtaining 
information on those memory kernels through MD calcu- 
lations of some realistic model of the supercooled liquid 
under consideration. Because the memory kernels are 
g — > limits of correlation functions of specific variables, 
they can, in principle, be computed from these MDcal- 
culations. Yet, these kernels, Eqs. ( |l^ , |2^ , p3| ) and (28a), 
are written in terms of the reduced operator R'{t). As 
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the latter has no easy formulation, this apparently re- 
duces drastically the value of the preceeding remark. We 
show, in the present Section, that the microscopic ap- 
proach of Section II allows for the determination of ex- 
pressions of rjbit), and i]T(t), the bulk and transverse vis- 
cosities, which coincide with the usual Green-Kubo for- 
mulae: they can be directly determined as the correlation 
functions of variables accessible in a MD calculation. 

Conversely, no such direct determination of F' (t) and 
fi{t) is possible; their Laplace transforms, can be ob- 
tained through the computation of the Laplace trans- 
forms of the correlation functions of other dynamical 
variables but the determination of r'(u;) and //(w) will 
be rather indirect, as we shall see; for technical reasons, 
we shall start with this second aspect and will turn, later 
on, to the determination of the viscosity kernels. 

Let us start by recalling that, in Eq. (|2^), the 'noise 
term' is equal to i?'(i)Qjj(q) so that one can rewrite this 
equation as: 

R'{t)Qo{ci) = Qo{ci,t)+ujlQo{q,t) 

+ I r'(t- s)Qo(q,s)rfs 



with: 



ro(q) = [2t,M) - r.M) - ryy{ci)]Nl2 . (41b) 



-A' / /^(t - s)ro(q, s)ds , (41a) 

^0 



When computing ^{t) or T'{t) through Eqs. 
t he q limit is taken, and ro(q) is 0{q) (see Eq. 
(24c)); the last term of Eq. (41a) may thus be dropped. 
Multiplying both sides of Eq. (41a), from the left, by Qo, 
performing a thermal average and a Laplace transform 
yields: 

-n'T'iu) - LT[{Q^\R'{t)\Q^)]{u:) 

= [ujI + LoT'iiu) - uj^)LT[{Q^{t)\Q^)]{uj) 
+ [lo~T'{u)]{Q^\Qo). (42a) 

As (QnlQn) = -(QolQo) = the n'^V'{uj) drops out 

of Eq. (42a): T'{lo) is not directly determined by con- 
sidering the 'noise term': its indirect determination is 
nevertheless possible through Eq. (42a) as: 



LT[{Q^{t)\Q)]{Lo) 



D{^) 



(42b) 



The l.h.s. of Eq. (421:) can be obtained from the cor- 
relation of Qo(i) with Qq. Nevertheless, it is simpler to 
write: 

LT[{Qo{t)\Qo)]{Lo) = LT[{Qo{t)\Qo)]{uj) 

= -uj''LT[{Q^{t)\Qo)]{Lo) + uj{Qo\Qo) . (42c) 



Eqs. (42b 42c) can be recast into the form: 



LT[iQo{t)\Qo)]icj) = — 

UJ 



1 - 



The fancy technique we have just used simply recovers 
Eq. (34) of Part I which was directly obtained from 
the phenomenological equations of motion. The latter 
have been microscopically derived in Section II, and this 
derivation implied the neglect of t he n oise term term 
R'{t). The above given proof of Eq. ( 42d ) can be consid- 
ered as a consistency check for the use of the 'noise term' 
to derive valuable results, a technique we shall now use 
to derive useful expressions for i]b{t) and ?7t(0- 

Before doi ng that, let us multiply the q ^ hmit 
of Eq. (41a), on the left, by 7r2o. Performing similar 
manipulations as above, one obtains: 

-n^fi{Lj) = LT[{Tr2o\R'{t)\Qo)]{u;) 

= D{uj)LT[{Qoit)\TT2o)Kio) 

~[LO-T'iio)]{n2o\Qo), (43a) 



where (7r2o|Qo) is equal to zero, due to Eqs. (21a) and 
(i). This yields: 



m(^) = - 

or, equivalently: 



LT[{Qo{t)\TT2o)]{uj) 



r(c.)=-^Lr[(Qo(i)k2o)](^). 



(43b) 



(43c) 



It should be noted that r{uj)/uj is of order 0{uj^^) for 
frequencies a; 3> wq; this is consistent with the 'sum rule' 
associated with {Qq\tt2o) = 0- f'i'^) is the function which 
couples the (longitudinal and transverse) phonon propa- 
gator to the light scattering mechanism via the orienta- 
tional part of t hese excitations, see Eqs. (31a) and (34a). 
Equation (43c) shows that r{t) can be directly obtained 
as the time derivative of the correlation function of Qoit) 
with 1T20, but that ii(t) is not directly accessible; it can 
be obtained only once r{uj) and D{cu)~^ have been deter- 
mined by the MD calculation. 



B. Expressions of rjT{t) and r]b{t) as time correlation 
functions 

Let us now use the same 'noise term' technique to 
express riT{t) and r]h{t) as auto-correlation functions of 
some dynamical variables. The 'noise term' of Eq. ( |25| ) 
is equal to i?'(f)7r2o(q), and, in the same q ^ limit, this 
equation simplifies into: 

R'{t)7T2o = 7r2o(i) - / Ai(i - s)Qo(s)ds . (44a) 
Jo 

Multiplying this equation from the left by tt2o, and per- 
forming the same manipulations as before yields, with 
the help of Eq. 1^ 



(42d) 



77s(w) = ir[(7r2o(i)k2o)](^) 
n 

-^i{uj){LoLT[{Qa{t)\^2o)]{oj) - (gok2o)}- 

(44b) 
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Using Eqs. ([43b|), (gla]) and (|), Eq. ( |44bD transforms 
into: 



iT[(^2o(i)k2o)](c^) + 



~ D{uj) 



(44c) 



From the definition of ririt), Eq. (31c), this equation 
reads: 



VT{t) 



knT 



(7r2o(i)k2o) 



(44d) 



Equation (44d) is the link between the usual Navier- 
Stokes approach to the dynamics of supercooled liquids 
and the most sophisticated approach of the present se- 
ries of papers, which takes explicitly into account the 
rotational motion of the molecules (the usual approach 
is rec overe d by formally putting ^{t) = 0). rjT^uj) (see 
Eq. ( ^4h| )) is the memory function which governs the 
transverse phonon pro pagator: within the Green-Kubo 
formalism of Eq. ( [44cj ) , this transverse viscosity is pro- 
portional to the correlation of the traceless part of the 
stress tensor, 7r2o, independently of the existence of a 
rotation-translation coupling. In other words, the pure 
center-of-mass viscosity, r]s{t), is not the quantity di- 
rectly measured by the correlation function of 1:20 ■ fjsit) 
must be deduced from the simultaneous determination 
of r]Tii-^j), and D{uj), quantities which can all be ob- 
tained, at least in principle, as correlation functions of 
some properly chosen variables, as we have just shown. 
Equation (44d) also proves directly, see Eq. (39c), that 



Imr]T{^) is always positive, being the Fourier transform 
of a auto-correlation function. 

The same type of technique can be used to determine 
r]b{t). In the q ^ limit, Eq. (|20b|) reads: 



R'it)p ^ Pit) - c^p{t) - R{t)[p - p{p\p)/{p\p)] , (45a) 
once Eqs. (|[) and (17c) have been taken into account. 



R\t)p = R{t)QnP. 



(45b) 



Equation (45b) introduces the variable QnP, which is the 
part of the pressure which is orthogonal to the density. 
Because of the existence of a Q operator on the left hand 
side of R'(t) (see Appendix B) which projects out the p 
variable: 



{QnP\Bl{t) = {p\R'{t) . 
Thus, Eq. ( pj| ) can be written as: 



ksT 



{QnP{t)\QnP) ■ 



(45c) 



(45d) 



Equation (45d) is the analog of the Green-Kubo formula- 
tion of the bulk viscosity within the usual Navier-Stokes 
formalism: as the rotation-translation coupling does not 
play a role in the bulk viscosity, this usual formulation re- 
mains exact in the more sophisticated present approach. 



Here a comment is in order. Since we did not deal with 
energy fluctuations in the projector, the correlation func- 
tion 77f)(t) decays to a non-zero constant even for times 
much larger than the structural relaxation time. It is 
therefore convenient to define a new correlation function 
fjbit) that vanishes at long times by an appropriate sub- 
traction. One can work out the constant from thermo- 
dynamic considerations and find: 



Vbit) = Vbit) + mn{c - c ) 



(46) 



where c is the adiabatic sound velocity. For the Laplace 
transforms this implies the relation: 



UJr]b{iLj) = ijjf}i,{Lo) + mn{c^ — c?) . 



(47) 



In all correlation functions considered so far, the bulk vis- 
cosity appeared only via the longitudinal phonon prop- 
agator. Using the preceeding equation PL{q,Lu) reads: 



~2 2 

c q 



f^[i^,(^) + LTiL0)], (48) 

mn 6 



which shows that the adiabatic sound velocity governs 
the propagation of longitudinal phonons. To simplify no- 
tations, in the remaining part of the paper, we shall drop 
the tilde again and treat c as the adiabatic sound velocity 
and ?7b(t) as decaying to zero for long times. 



V. COMPARISON WITH PREVIOUS 
THEORETICAL APPROACHES 

A. Introduction 

The discussions performed in and |^ have made 
clear that the set of Eqs. ( p6| ) and ( p9| ) are convenient 
tools to describe the light scattering spectra of molecular 
liquids, in their normal and in their supercooled states, 
when those equations are supplemented by the dielectric 
model of Eq. (^0|). Indeed as soon as the four memory 
functions riij{t),i]s(t), p{t) and T'{t) are mimicked by rea- 
sonably decreasing functions (characterised, inter alia, by 
relaxation times, r, that increase with decreasing temper- 
ature) the most characteristic features of the VH spectra 
can be described: 

- The back scattering spectrum is mostly characterised 
by a broad high-frequency libration mode, in the vicinity 
of a frequency a;o/27r, and by a low- frequency central 
mode, the line width of which decreases upon cooling. 
Both features can be approximately reproduced by Eq. 
( |34d| ) with the help of the expression of D{lu) given by 
Eq. (|3^), as soon as a reasonable T'{uj) is chosen. 

- The shape of the q-dependent part of the VH spec- 
trum has been discussed in detail in ||l[. It was shown 
that Eqs. (|4^ and ( |34b| ) allowed to adequately describe 
the existence of a Rytov dip [l4| in a normal molecular 
liquid, this dip being a very narrow central peak, wave- 
number and scattering-angle dependent, which is sub- 
tracted from the much broader central mode. This dip 
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appears in the high-temperature regime when, for all the 
frequencies of the central mode, lot <C 1. The ujt ^ 1 
regime, which is characterised by the appearance of the 
Brillouin spectrum of a transverse propagative mode, is 
also well described by these equations, provided reason- 
ably decreasing functions are also taken for the three 
remaining memory functions. In particular, the trans- 
verse sound velocity, characterised by the plateau value 
of ujr]T{uj) at frequencies 1 ^ ujt ^ uqt is decreased by 
the coupling of the molecular orientation to the trans- 
verse phonon throu gh fi {uj)'^. Finally, in view of the form 
of 77l(w), see Eq. (31c), the same is true for the sound 
velocity of the longitudinal phonons. 

In this Section, we shall cornpare the results which 
can be obtained through Eqs. (|2|) and ( p9|) with those 
resulting from the two other papers (or series of papers) 
already mentioned which make use, in different ways, of a 
Mori-Zwanzig technique to describe the liquid dynamics. 



B. The Andersen- Pecora approach 

The Anderson and Pecora approach Q was only used 
to study the VH spectrum of a molecular liquid at high 
temperature |21 . Indeed, their work was devoted to 
the explanation of the Rytov dip IQ; and their anal- 



ysis made use of a dielectric fluctuation model identical 
to the one of the present paper. 

In the work of Andersen and Pecora, the mass den- 
sity, the mass current and a second-rank tensor propor- 
tional to Qij were the sole 'slow variables' of the theory, 
within the usual Zwanzig-Mori distinction between 'slow' 
and 'fast' variables. In other words, they implicitly as- 
sumed that Qij had a much faster dynamics than Qij, 
so that the former could be treated on the same footing 
as the other fast variables. Furthermore, they performed 
a Markov approximation on all the retarded interactions 
that needed to be taken into account, which is equivalent 
to taking the ujt ^ 1 limit of the corresponding kernels. 

A summary of the result of their theory, within this 
Markov approximation, can be found in the book of 
Berne and Pecora |l^. The corresponding equations 
read, with notations adapted to the present paper: 

Try (q, t) = -r'/iTy (q, t) - iT'.^Q,, (q, t) , (49) 



Q^3 (q, t) = -iP^iTy (q, t) - T22Q^J (q, t) . (50) 



Here, the kinetic coefficients r"2,r22 are real quantities, 
whereas r']^2ir2i are purely imaginary, and are related 
by Onsager's principle. Equation ( [49| ) makes it clear 
that 7ry (q, i) depends, here, linearly on Qij{c[,t), and 
not on its time derivative, as is the case in Eq. (25), 
while Eq. (^0|) does not contain a second time derivative 
of Qy (q, t), contrary to Eq. ^ 



The form of the Zwanzig-Mori technique, see Eq. (14) 



variables of the problem (time dependent generalisations 
of T'(i , T'i2 ^iid 1^22 ) , in terms of a reduced time evolution 
operator. Calculating from the corresponding equations 
of motion the VH spectrum and comparing the results 
with Eqs. (34d, 34e), one can express the three Andersen- 
Pecora kernels as functions of ris{uj), ^{u!),r' {oj), A' and 
loq. One can thus study their lot <C 1 and lot 3> 1 
regimes. This study will show that the lot <C 1 limit 
gives reasonable results, which are, as expected, in line 
with the Andersen-Pecora Markov approximation. Con- 
versely, the LOT 3> 1 limit yields a complicated behaviour 
for the same three kernels which cannot be easily mod- 
eled. This will make the Andersen-Pecora method inap- 
propriate for the study of a molecular supercooled liquid, 
as we shall now see. 

Indeed, using the same technique as in Section II, one 
can derive the equations of motion related to H.y (q, t) 
and Qij (q, t) when one restricts the variables to the 
Anderson-Pecora set. This means that, e.g. the pro- 
jection operator P of Eq. dl^ ) has been replaced by: 



P = \QM)^{QkdcD\ + \p{ci))^ip{q)\ 

jLitJ 111 D 

+ |Jfc(q))^(Jfe(q)|. (51) 



P leads to the new orthogonal projector Q — 1 — P, in 
terms of which a new reduced time evolution operator 
R{t) can be defin ed. One then finds that the equation 
for p(q, t), Eq. ( ^Ob| ) , is not modified, except for the 
change of R'{t) into R'{t), which formally changes ?7b(t) 
into fjh (t) . Evaluating the corresponding Green-Kubo re- 
lation reveals that 77b (t) = rib{t). Since the pressure fluc- 
tuations are irrelevant for the VH spectrum discussed in 
0, we do not discuss further those aspects. Conversely, 
the equation for TTy (q, t) now turns out to be: 

T^tji<i,t) = - / rjsit- s)Tij{q,s)ds 
Jo 



+ / X{t — s)Qij{q, s)ds + noise , (52) 



f/sit) being the complete analog of ris{t)- Equation (|52 
contains a linear term in Qij{q,t), as in Eq. (^), and 
not in Qij{(l,t), as was case for Eq. ( p5| ) while the cor- 
responding retarded interaction is expressed by: 



{Q;j\R'{t)\TTkl)-^ = -A(t)Ay-fc, 



(53) 



where Qij and ivki are respectively odd and even with 
respect to time inversion; thus X(t) is an odd function of 
t, contrary to all the memory functions considered up to 
now in the present paper. The most important change 
arises, nevertheless, from the fact that the equation of 
motion for Qij{q, t) has to be derived from: 



used in the present paper allows to derive precise expres- 
sions for the three relaxation kernels associated with the 



dtQij{q,t) = Qijiq,t) , 



(54) 
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an equation which replaces Eq. (12). A calculation simi- 
lar in every respect to the one performed below Eq. ( |2^ ) 
yields: 

A' /■* 

Qij{q_,t) = 2 / A(t - s)ry (q, s)(is 

'^o Jo 
t 

M{t — s)Qij (q, s)ds + noise , (55) 



with: 



(56) 



Equations (52) and ( |55|) are, obviously, the non- 
Mar kovian form of Eqs. ([49|) and (pO[). Ignoring temper- 
ature fluctuations, the previous relations are exact, and 
allow us to relate the memory kernels M{uj),fis{iv), X{iu) 
to the ones already used in this paper by deriving from 
Eqs. (52l5^), through the same methods as used in [Q, 



the expression of the VH spectrum and comparing it with 
Eq. (48) of Part I. For M{lo), this can be done by simply 
computing the correlation function of Q±±' which is re- 
sponsible for the pure back-scattering spectrum (see Part 
I for a definition of the geometry used) . Solving Eq. ( |55| ) 
in this simple case yields: 



ir[(Qii,(q,OIQ±±'(q))]M 



uj - M [uj) 



(57) 



Comparison with Eq. (34) of Part I, with the help of Eq. 
(|6a|), leads to the relation between T'{u}) and M{uj): 



M{uj) 



(58) 



In order to gain some insight into the Markov approx- 
imation, a priori valid at high temperatures, let us dis- 
cuss the properties of M{uj) upon cooling the system. To 
simplify the discussion, we consider a Maxwell model for 

r'(^): 



r (w) = 17 -f 



1 -I- iuiT 



(59) 



which mimics its frequency dependence for frequencies 
much lower than typical liquid frequencies, i.e. for 
UJ 't^ loq. Then, all the fast processes are hidden in a 
weakly temperature-dependent background, ij, whereas 
the temperature-sensitive structural relaxation is mod- 
eled by a decreasing exponential in the time domain, 
corresponding to a temperature insensitive amplitude, 
Fq, and a relaxation time, t, that increases by orders 
of magnitude upon supercooling the liquid. 

At high temperature, i.e. for lut <^ 1, under the pos- 
sible conditions uj, 7, ujq ^ F§t, which simply imply that 
even at high temperatures FqT is subst antia lly larger 
than the width of the central peak, Eq. (34d), one ob- 
tains: 



(60) 



This is the Andersen-Pecora result in its Markov approx- 
imation M{uj) = jF22 with F22 ~ T~^. Conversely, at 
low temperature, and under the same conditions, M{uj) 
will be approximated by: 



M{uj) 



lUJJ 



F2 
^ 



(61) 



at low frequencies, M{uj) is a real quantity proportional 
to UJ and independent of the relaxation time. The van- 
ishing of the imaginary part of M{uj) at low frequencies 
implies that the area of M{t) cancels at low tempera- 
tures. Whereas the shape of M{t) is clearly model de- 
pendent, the cancellation of areas of M{t) is a general 
feature of supercooled hquids. Contrary to T'{t), there is 
no step process in M(t) with a diverging time scale upon 
cooling. It is thus fruitless to try to model the tempo- 
ral evolution of M{t) since the common features of an 
increasing structural relaxation time are masked in this 
approach. 

We can similarly evaluate the Andersen-Pecora mem- 
ory kernels J7s(w), and X{uj) by solving the dynamics for 
the variable Q±\\{ci,t) that also contributes to the VH 
scattering. Using the methods of one obtains: 



XT[(Q^,[(q,<)|g^„(q))](c^) 
A' \{uj) 



M{uj) 



M{uj) 



2 22 



UJ — q'^rjT{uj)/mn 



with the transverse viscosity given by: 

A' A(w)2 



7?t(w) = ^.{UJ) 



M{uj) 



(62) 



(63) 



A comparison betw een the second term of the r.h.s of Eq. 
( |6^ ) and Eq. ( |34e| ) yields the relation between the two 
sets of memory kernels: 



X{uj) 



^oM(^) 
v~T'{uj) 



risiuj) = ris{uj) 



-ili{uj)M{uj) 



_A^ji(ujf_ 

cj-F'(w) 



(64) 



UJ, 



,^i{ujYM{uj) 



(65) 







The high-temperature limit of M{uj), Eq. (60), yields, 
with a Debye model for ^{uj) and rjs{uj), the Markov lim- 
its obtained in [Q: \{uj) becomes an imaginary number 
independent of r, while f}s{uj — s- 0) is the sum of two 
terms, both imaginary and proportional to the relaxation 
time. One originates from 77^(0; 0), while the second, 
negative, is the a; — > limit of K' ii{ujY M (uj) / uS^; this 
explains why, in the Andersen-Pecora approach, the vis- 
cosity, —irjs{uj 0), is the sum of two positive terms. 

Conversely, the difficulty of an a priori modeling of 
M{uj) transfers to the two other memory kernels, fis{uj) 
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and X(uj). This explains why a nai've modeling by simple, 
Debye-like, relaxation functions, consistent both with the 
high-temperature Markov results and the different time 
reversal symmetries of the kernels, is unable to yield cor- 
rect physical results in the supercooled regime. Appendix 
D shows, indeed, that the low-temperature limit of such 
an attempt leads to the existence of transverse propaga- 
tive modes coupled to molecular orientation motions, but 
this coupling increases the sound velocity instead of de- 
creasing it. 

Summarizing this part, we have shown that a Mori- 
Zwanzig procedure as applied by Andersen and Pecora 
Q allows to derive constitutive equations through which 
the light scattering problem can be properly formulated 
as long the frequency dependence is kept on a formal 
level. Conversely, when one expresses these memory ker- 
nels in terms of those obtained in Section II, one discovers 
that their modeling as time-dependent memory kernels is 
extremely difficult. This problem can be circumvented, 
if one models them directly, as inspired by Eqs. (58), 
(U), (|65|), but this procedure is equivalent to consider- 
ing T' (oj) , ii[lu) , and ?7s(w) as fundamental quantities. 



C. Comparison with the general expressions for 
light scattering 

The expressions for the VV and VH intensities ob- 
tained in P art I and discussed again in Section IV, Eqs. 
(31a) and (34a), have been obtained under the physical 
assumption that the fluctuations of the local dielectric 
tensor, 6eij{c[,t), could be expressed through Eq. (30), 
i.e. depend, in first order, only on the density and the 
orientational fluctuations. Conversely, the expressions 
obtained in 0] did not make use of a specific form for 



Seij (q, t) . We shall show, in this last part of Section V, 
that those two expressions are, indeed, a specialization of 
the general results obtained in ||] . We also discuss the re- 
spective merits of these two complementary approaches. 

The basic idea of was to express the finite wave- 
number fiuctuations of the dielectric tensor in terms of 
the long-wavelength limit of the two special linear com- 
binations: 

SOo(q) = [exxicO+eyyicO+ezziq)]/^, 

i2o(q) = [2£zz{q) - exx{<i) - '^vv{<l)]/VT2 ■ (66) 

Within the present light-scattering model, these two 
quantities reduce to the contributions of the density and 
the orientation fiuctuations, respectively: 



that have been adapted to the current paper, as: 

4_, m'^v'^ / dspo \ ^ 



9soo 



Loavv{^) 



dp / T 

e(c.)2LT[(e(q,t)|e(q))](c.) 
f dsoa\ 

T 



LTMq^tMqMco) 
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m V 



LT[(p(q,t)|e(q))](^) ^68) 



Similarly, the VH intensity was expressed as: 



/y//(q,w) ^ Im< T{uj) 



-—cos^ -[ujavH{uj)]^'m^v^PTiq,uj) } . (69) 



The transverse phonon propagator is given directly in 
terms of the transverse viscosity, 77^(0;), as in Eq. (^4^). 
As already alluded in the Introduction, the price to be 
paid for these two general results was the introduction of 
ten frequency-dependent quantities, namely, the scalar 
background spectrum S{lu), the tensor background spec- 
trum T{uj), the two Pockels' coupling functions, ayvi'^) 
for polarised and av^r(w) for depolarised scattering, and 
the temperature coupling ^(w), while the three hydrody- 
namic correlation functions related to the density, p(q), 
and kinetic temperature, 0(q), were expressed in terms 
of the transverse viscosity, 77^(0;), the longitudinal vis- 
cosity, •qL{uj), the heat conductivity, \{lo), the dynamic 
specific heat, cv{lo), and the tension coefficient, f3{uj), 
respectively. The singular hydrodynamic behavior man- 
ifested itself explicitly in the three correlation functions 
just mentioned. 

Let us relate those quantities to the one derived in 
the present paper and demonstrate numerous simplifica- 
tions that occur in the density-and-orientational-decay- 
channels-only model. First, the temperature coupling is 
given by: 



cv{l^) cviuj) LT[iQsoo{t)\QeP)]{u;) 



cv 



knT^ 



, (70) 



soo(q) = ap(q) , 



(67a) 



where ^ is a linear function of the energy fluctuations in 
the g — > limit and Q = 1 — P is a projection operator 
orthogonal to the standard Kadanoff-Martin projector, 
P. The latter projects, in ||], on five variables, p(q) and 
J(q), as in the present paper, and on the temperature 
fluctuations, T(q), not introduced here, and proportional 
to the energy fluctuations, e(q). 



t2o(q) = &Qo(q)- (67b) 
In H, the VV spectrum was expressed, using notations 



Hr(q)) 



cy 



(T(q)|+0(92), 



(71) 
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Because energy fluctuations are not considered in the 
present approach, ^ = 0. Also, Soo(q) is proportional 



to p(q), Eq. ( |67ci| ) , so that Qsqq — 0,i^(w) = 0. 

Second, in ||5|, the scalar and the tensor background 
spectra are defined as: 



Situ) 



-LT[(sooWIQsoo)](^), 



Tiio) = iT[(t2o(i)|i2o)](w). 
The preceeding results imply: 

Situ) = 0, 

TH = b'LT[{Qo{t)\Q„)]{iu). 



(72a) 



(72b) 



(73a) 



(73b) 



Third, from Eq. ( |67a| ), {dsoo/dp)T — a while the dy- 
namic Pockels' coupling functions are given, in ||^, as: 



avHH = LT[{Tr2oit)\t2o)]{Lo)/m'v\ (74a) 

ayv(w) = -avH[^^) H 

3 UJCv(i^) wcy 

-LT[{p{t)\Qsoo)]{Lo)/mW (74b) 
As ^(oj),^ and Qsqq are all equal to zero, one obtains: 



which is here a simpl e definit i on of uj^ I Djuf) . In the same 
manner, from Eqs. (74a, 67b ,43c, 28b| ), one obtains: 



mn 



(78) 



where, similarly to the case of luq/ D{LL!),r{uj) is simply 
defined through the Laplace transform of (7r2o(i)|i2o), Eq. 
(74a). One sees that Eq. ( [76|) has been cast into a form 
identical to En. (31a), while a si milar identification holds 
between Eq. @) and Eq (|34a| ). 

The reduction of Eqs. (|68|) and to Eqs. ( |3Ta| ) and 
4a^) shows the comparative interests of the approach of 
and of the present one. The method of |^ does not 
depend on the system under study and allows for tem- 
perature (or energy) fluctuations: as soon as the scat- 
tering model, Eq. (^0[), is introduced, and the energy 
fluctuations are neglected, the equations of |5| reduce to 
those of the present model, depending on four functions 
T{uj), avni^), VL{i^) and riT{uj), which are undetermined 
in this framework. Conversely, the more restricted ap- 
proach developed in the present series of papers gives 
precise definitions of these four quantities in terms of 
more fundamental memory kernels T'(lu), /i(w), ilbi^) and 
77s (cj), and also gives the relationships through which the 
four first functions are related to the second ones via the 
constants A' and luq for which definitions can be given. 
Yet the restricted approach has its own price to be paid: 
one has to start the whole work again if additional vari- 
ables need to be introduced into the model. 



(74c) 



Let us look at the results for the VV light scattering 
intensities. One observes that terms involving dynamic 
correlation functions of the kinetic temperature evaluate 
to zero. One is thus left with the density correlation 
functions which, if energy fluctuations are ignored, reads 
in agreement with Eq. (29a) of Part I: 



ir[(p(q,i)|p(q))](c.) = 



uj^ — [c^ + UJ7]L {uj)/mn\ 



(75) 



Equation ( [75| ) allows to group the terms proportional to 
rri^v^ and one ends up with: 



Ivv{<l,uj) = Imi ■^T^i'^) 



a - -ujavH{to) 



^2 _ q2^2 q'^ujr]L{i^)/mn 



Tjuj) and avni^) being deflned through Eqs. ( [73b| ) and 
(74a), respectively. Furthermore in agreement with Eq. 
(42c) and Eq. (34) of Part I, T{uj) can be written as 



D{uj) 



(77) 



VI. SUMMARY AND FINAL REMARKS 

In a liquid formed of rigid molecules, the dynamics of 
the system has to take into account both the motion of 
the molecular centers of mass and the orientational mo- 
tion of the molecules. In the long-wavelength limit, the 
flrst one gives rise to the hydrodynamic modes, to which 
the orientational motions are partly coupled, while this 
orientational dynamics also gives rise to motions that are 
wave- vector independent in the same limit. ^ ^ and |^ 
proposed a phenomenological set of equations to describe 
this coupled dynamics in the case of linear molecules, and 
a phenomenological expression for the local fluctuation 
of the dielectric tensor: this fluctuation was expressed in 
terms of the density and orientational variables entering 
the dynamical equations. 

The original objective of the present paper was twofold 
- one was to provide a complete, microscopic, derivation 
of these dynamical equations; the second was to compare 
the expression for the light scattering intensities resulting 
from these equations with those obtained with two other 
approaches 

Both goals have been achieved. On the one hand, the 
use of a Zwanzig-Mori formalism has allowed to com- 
pletely derive these dynamical equations; in the course of 
this derivation, we have obtained the microscopic expres- 
sions of the two parameters and of the four memory func- 
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tions entering those equations. On the other hand, the 
comparison with the two other Zwanzig-Mori approaches 
has also brought important results. One of them is re- 
lated to the choice of Andersen and Pecora jj] of not 
including Qy (q) in their set of variables. This choice, 
which is sufficient at high temperature, when the Markov 
approximation can be made on the corresponding mem- 
ory kernels, turns out to be inappropriate for the study 
of supercooled liquids: at low temperatures (3ij(q, t) is 
as 'slow' a variable as Qij{q,t). The most important 
consequence of the absence of Qij (q, t) in the set of se- 
lected variables is the change in the equation of motion of 
Qij (q, t) : it transforms it from a second order differential 
equation with a memory kernel acting on Qij(q) into a 
first order differential equation with a memory kernel act- 
ing on (5ij(q). The formal neglect of (q) in the set of 
'slow' variables is possible but the corresponding memory 
kernels have a non-trivial time evolution, which cannot 
be predicted without using the results of the present the- 
ory. The second result is that the present formulation of 
the theory is, indeed, a reduction of the general theory of 
which can be derived from simplifications consistent 
with the phenomenological model of the dielectric tensor, 
and with the restricted set of variables used here. 

The present Zwanzig-Mori approach also led to two im- 
portant bypr oducts. One is the existence of conditions, 
Eqs. (^,38b), which have to be fulfilled by the Laplace 
transforms of the memory functions. These conditions 
are important because they are sufficient to insure that 
all the light scattering intensities will be positive, what- 
ever the frequency, within the phenomenological model 
of the fluctuations of the dielectric tensor used here. A 
second byproduct is the set of Green-Kubo formulae we 
have derived in Section IV-C: we have shown that the 
correlation functions of some variables, not experimen- 
tally accessible by light scattering techniques, but which 
may be numerically obtained from MD computations of 
models of these molecular (supercooled) liquids, give ac- 
cess to definite combinations of the Laplace transforms 
of these memory functions. This is a possible way of 
obtaining an information on them. 

Some results of the present paper provide a direct help 
to the experimentalists, when analysing the light scat- 
tering spectra of molecular supercooled liquids formed of 
rigid linear molecules, or of molecular liquids for which 
such an approximations is reasonable. One of them is 
the already mentioned necessary conditions on the mem- 
ory functions. A second is that these functions exhibit 
the characteristic features of structural relaxation, e.g. 
rapidly increasing relaxation times upon lowering the 
temperature. Yet, the functional form of these mem- 
ory kernels remains undetermined within this framework, 
except for the conventional analytic properties. Usually 
it is not possible to directly extract the frequency de- 
pendence of of the memory kernels from light-scattering 
experiments. Rather one has to rely on empirical func- 
tions and adjust a small number of parameters to ob- 
tain a reasonable description of experimental data. As 



a further step, one can supplement these empirical func- 
tions with features inspired from theoretical considera- 
tions, e.g. the fast /3-process as discussed in the context 
of mode-coupling theory 1^. 

Time resolved optical spectroscopy of the same molec- 
ular liquids has recently developed into an important 
tool; this is particularly the case for the impulsive stimu- 
lated thermal scattering technique (ISTS) mentioned in 
the Introduction 0, |^. The most important part of 
the new information obtained from these measurements 
is derived from the coupling of the heat diffusion pro- 
cess with the stimulated hydrodynamics mode. We have 
not incorporated, in the microscopic derivation of the 
dynamical equations, a local temperature as a pertinent 
variable, contrary to what has been done in ||^. In order 
to properly exploit the information contained in these 
ISTS experiments, the whole procedure developped in 
the present paper has to be repeated with the inclusion 
of the variable(s) describing the local temperature of the 
supercooled liquid. It has to be found if this general- 
isation will require more memory functions than could 
be anticipated from a phenomenological extension of the 
full set of Navier-Stokes equations to the case of a super- 
cooled (memory function aspect) molecular (inclusion of 
the rotation-translation coupling and of the molecular 
orientation dynamics) liquid [l5[| . 
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APPENDIX A: STATIC AVERAGES 

The Hamilton function of identical, interacting, sym- 
metric tops reads; 

H = J2T^ + V{{R^,(j)^,e^}), (Al) 

a 

where the kinetic energy of the a-th molecule is given by: 



p2 

Ta = — 

2m 



2/sin2 Qa 



^ + ^. (A2) 



21' 



Here R^, 0ct, Oq, denote the center-of mass position 
and the Euler angles of the molecule following the def- 
inition of ||l^, and Pa,Pa4>,Pae,Pati! the corresponding 
canonical momenta. The moments of inertia are denoted 
by /, /' for rotation perpendicular to and around the 
molecule's axis of symmetry. The potential energy of 
the interacting molecules is denoted by V. Note that, 
due to the symmetry, the interaction does not depend on 
the Euler angles ipa. 
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The orientational current, Qy(q) = iCQij{q) 
{H,Qij{q)}, then sphts naturaUy mto two parts: 



N 



,(q) = iV-V2^ 



a=l 
N 



(A3a) 



where the first term corresponds to the translational mo- 
tion of the center of mass and the second term describes 
molecular reorientations. 
By definition: 



(^(q)|s(q)) 



dTe-"^''''^6A{q)*dB{q) ,(A3b) 



where dT — YladIlad(l)adQadiJ-'adPadpaci)dpaedpa^ de- 
notes the canonical phase space volume element. Let us 
compute {Qij{q)\Jk{q))- In this t herm al average appear 
two types of integrals (see Eq. ( A3a )): one involving 
(q • Pa)P/3k, and the second (iLuaiUpj)Ppk- As 



(A3c) 



One checks that the Jacobian is sin0Q,, while the corre- 
sponding part of the kinetic energy reads — {J^x + 
)/2J + . Then the partition sum is given by 



Z = 



j Y{d'Ra,d(j)adcosQadipadPad3aer"'^'''^ ■ (A6) 



Thus, averaging over the angular momenta is just gaus- 
sian and averaging over the Euler angles amounts to aver- 
aging over the usual Haar measure of the rotation group. 

In order to calculate the long-wavelength limit of the 
auto-correlation function of the orientational currents, 
it is sufficient to calculate it for one of its components, 
say Qzz{q)- Since = cos9q and 8q = {H,Qa\ = 
PaG)/!, one finds, for q — > 0: 



N 



Qzz{q) = -jN-^/'^^y^x^uiiP^ 



a=l 

+ Jay COS V'a] COS Qa siu Qc 



(A7) 



Then the long-wavelength correlation function of the 
orientational current can be evaluated: 



involves only the angular variables, the only part depen- 
dent on t he li near mo mentum in the second terms reads 
(c.f. Eqs. ^and^): 

e-^V^'-'^^^P,,dP,, 



0. 



(A3d) 



One is thus left with the contributions of the first term. 
They read: 



N 
a,/3=l 



iq • Pa 



m 

N 



Q = l 

N 

= -iqkkBTN^^'^{{uaiUaj 



3 ^ 



(Q..(q = 0)|g..(q = 0)) = 
4 



P 



^ '^{\jax Sin?/>Q -I- Jay 



COS t/iq 



cos <da sin Oc 



cos 9„ sin Qr 
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(A8) 



a=l 



(A4) 



Comparison with Eq. (|^) yields for the ideal gas libration 
frequency fl'^ = 2kBT/5I. 



APPENDIX B: OPERATOR IDENTITY 

The time evolution operator, R{t) = exjp{iCt), may be 
split into two parts R{t) = Rp{t) + Rqit) with Rp{t) = 
R{t)P, RQ{t) = R{t)Q. From the equation of motion, 
dtR{t) — R{t)i£, one finds: 



dtRgit) = Rp{t)iCQ + RQ{t)iCQ . 



(Bl) 



where the last but one equality originates from averaging 
over the gaussian variable, Pq, and the last one from the 
rotational symmetry of the problem. Equation (A4) has 
been reported as Eq. (H) in the body of the present paper. 

The kinetic energy expressed in terms of canonical mo- 
menta depends explicitly on the Euler angles, hence the 
evaluation of thermal averages is quite involved. This 
can be avoided by eliminating the canonical momenta in 
favour of the angular momenta p^ : 

Pa(t> = —J ax sin Qa COS Tpa + Jay sin 8q siu Ipa 

+ Jaz COsQa , 
Pae = J ax sin -0Q -I- Jay COS Ipa , 

Palp — Jaz ■ (-^5) 



The solution of Eq. (Bl) can be expressed in terms of 
Rp{t) as: 



RQit) 



iCQt 



Rp{s)iCQe'^^'^*-'Us. (B2) 



Furthermore, because R'{t) — Qe incorporates the 
projection operator Q, one easily finds, by e.g. expan- 
sion of the exponential, that R'{t) may be written in the 
symmetric form: 



R'{t) = Qe"^^'^'Q. 
Collecting terms, one arrives at Eq. (p^ 



(B3) 
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APPENDIX C: PROPERTIES OF THE MEMORY 
FUNCTIONS 

The memory kernels of the type (v4|i?'(i)|A) exhibit 
the same mathematical properties as auto-correlation 
functions, viz. for complex frequencies in the lower 
half plane, their Laplace transform is analytic with non- 
negative imaginary part. A non-rigorous proof can be 
adapted from Berne and Pecora ||T^. Since R'{t) = 
Q exp(iQCQt)Q, we formally introduce a complete set 
of eigenfunctions of the hermitian (with respect to the 
Kubo scalar product) operator Q£Q: 



APPENDIX D: TRANSVERSE PHONONS AND 
THE ANDERSEN-PECORA APPROACH IN THE 
LOW TEMPERATURE LIMIT 

If one makes the (incorrect) supposition that the mem- 
ory kernels of the Andersen- Pecora approach can be mod- 
eled by Debye relaxation processes consistent with their 
high-temperature Markov approximation and their time 
reversal symmetry, this yields: 



M(uj) 



1 



T 1 + iuJT 



(Dl) 



QCQ(f)\ = A(/)A , 
where all eigenvalues are real. Thus we can write: 



(CI) 



{A\R'{t)\A) = {QA\e'^'^^'\QA) 

= ^(Q^|0A)e*'*(0A|Q^). (C2) 

A 

The Laplace transform yields: 

LT[{A\R'{t)\A)]{u:) = ^ ^|(QA|0,)|2 , (C3) 

with complex frequencies in the lower half-plane. Since 
all the poles are located on the real axis, the Laplace 
transform is analytic for uj — ^ — ie,e > Q. Furthermore: 



ImLT[{A\R' {t)\A)]{u) 



E 



{a - A)2 + C 



■|(QA|0a)P>O. (C4) 



In particular, provided the limit e \ exists, one obtains 
for real w. 

ImLT[{A\R'{t)\A)]{u) = ^7r^(c^ - A)|(QA|0a)P > 

(C5) 

Consider now a collection of phase space variables 
Ai,i — 1 , . . , Z of identical time inversion parity. Then 
the real symmetric matrix ImLT[{Ai\R' {t)\Ajy\{uj) is 
positive semi-definite: since for arbitrary real numbers 
yi,i = 1,..,Z the spectrum of the autocorrelation func- 
tion oiY = X]i=i Di-^i is non-negative, one finds: 



^ V,V,ImLT[{A,\R'{t)\A,)]{uj) > 0, (C6) 

which implies the property. For frequencies w — > 0, one 
obtains Onsager's relations, viz. the matrix of the kinetic 
coefficients i s sy mmetric with non-negative eigenvalues. 
Hence, Eq. (C6) can be interpreted as the proper gener- 
alisation of Onsager's relations to finite frequencies. 



A(^) 



ifl + /2 



1 -I- iuJT 



(D2) 



1 + iujT 



/i,/2>0. (D3) 



The special form proposed for Eq. (D3) derives from the 
fact that A(i) is an odd function of time. If we suppose 
it to be the time derivative of /(t), the auto-correlation 
function of some variable, /i is its f = value and we 
have chosen for its late-time evolution a smooth Debye- 
like behaviour. 

Let us admit that, in the cut 3> 1 regime, the value of 
A is smaller than l/v^. Equation (Dl) then yields a VH 
backscattering spectrum with a pscudo-Lorentzian line 
shape and a line width approximately equal to r~"'^A2(l — 
2^2)"^/^. For uj larger than this line width, one can write 
for the g-dependent part of Eq. (|62|) : 



I{q,uj) = ~q 

xlm 



A' 



2 9 9 

m V 



1 



Re 



Ll!X{lu) 

' - M{lu) 



In this UJT ^ 1 limit, this reads: 

A' 



(D4) 



2 9 9 

m V 



UJ 



(/l - hf 



xlm 



1 



— q^cl^ 



(D5) 



e being a small positive quantity. Equation (D5) does 
represent the Brillouin spectrum of a transverse phonon, 
but the square of its the velocity is given by: 



Crp 



mn 



-^{fl-f2? 

mnujQ 



(D6) 



In the spirit of the Andersen- Pecora approach, fisit) 
represents the contribution of the molecular center-of- 
mass motion to the shear viscosity (see, nevertheless, 
the remark below Eq. (|65|)). Within the same spirit, 
[?7j/mn]^/2 should represent the contribution of the same 
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motion to the transverse sound velocity. As announced, coupling, but this renormalisation leads to an unphysical 
the bare sound velocity, [^°/mn]^/^, is renormalized by a increase of the sound velocity, instead of the expected 
term in A'/wq, the signature of the rotation-translation physical decrease. 
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